Carbonates in space — 
The challenge of low temperature data 



i< ; Th. Posch, A. Baier 

Institut fur Astronomie, Tiirkenschanzstrafie 17, A-1180 Wien, Austria 
CN . posch@astro .univie . ac . at 

and 



OO 



(N 
> 

On 

m 



X 

3 



H. Mutschke 



Astrophysikalisches Institut, Schillergasschen 2-3, D-07745 Jena, Germany 

^! and 
O 

c3 '. Th. Henning 



Max-Planck-Institut fur Astronomie (MPIA), Konigstuhl 17, D-691 17 Heidelberg, Germany 



ABSTRACT 



Carbonates have repeatedly been discussed as possible carriers of Stardust 
emission bands. However, the band assignments proposed so far were mainly 
based on room temperature powder transmission spectra of the respective min- 
erals. Since very cold calcite grains have been claimed to be present in protostars 
and in Planetary Nebulae such as NGC 6302, the changes of their dielectric 
functions at low temperatures are relevant from an astronomical point of view. 

We have derived the IR optical constants of calcite and dolomite from re- 
flectance spectra - measured at 300, 200, 100 and 10 K - and calculated small 
particle spectra for different grain shapes, with the following results: i) The ab- 
sorption efficiency factors both of calcite and dolomite are extremely dependent 
on the particle shapes. This is due to the high peak values of the optical con- 
stants of CaCOs and CaMg[C03]2- ii) The far infrared properties of calcite and 
dolomite depend also very significantly on the temperature. Below 200 K, a pro- 
nounced sharpening and increase in the band strengths of the FIR resonances 
occurs, iii) In view of the intrinsic strength and sharpening of the ~44 /im band 
of calcite at 200-100 K, the absence of this band - inferred from Infrared Space 
Observatory data - in PNe requires dust temperatures below 45 K. iv) Calcite 
grains at such low temperatures can account for the '92' /im band, while our data 
rule out dolomite as the carrier of the 6^-65 /mi band. 

The optical constants here presented are publicly available in the electronic 
database http: //www. astro .uni-jena.de/Laboratory/OCDB. 



Subject headings: circumstellar matter — infrared: stars — methods: laboratory - 



1. Introduction 

In the 1970s, carbonates have been pro- 
posed as carriers of mid- infrared features 
detected in astronomical spectra, such as 
the 11.3 /mi emission band seen in NGC 
7027 and in the 'Red Rectangle' (Gillet, 
Forre st fc Merr i ll 197 3; Bregman & Rank 
1975: iBregmlr] (jl977h ). 

In parallel to that, reflectance spec- 
tra of calcite (CaCOa) and magnesite 
(MgC O.s) have been measure d in the labo- 
ratory. Hellwege et al. ( 1970l ) did the mea- 
surements in the whole mid- and far-IR 
range for calcite, mag nesite an d dolo mite 
(CaMg[C0 3 ] 2 )), while |Penman| (jl97fih ex- 
amined calcite and magnesite in the 5- 
50 region. All these measurements 
have been carried out at room tempera- 
ture. 

Carbonate grains have been fou nd in in- 
terpla netar y dust particles (IDP s ). ISandford 
(119861 ) and llbmeoka fc Buseckl (119861 ) re- 
ported the identification of carbonate min- 
erals in IDPs both by a chemical dis- 
solution experiment and by transmis- 
sion electron microscopy . As shown by 
Tomeoka fc Buseckl (119861 ) , the mineralogy 
of the IDP carbonates is essentially that 
of breunnerite - an Mg-rich solid solution 
of magnesite and siderite (FeCOs). Solid 
solutions of Ca- and Mg-carbonates are 
rare in the IDPs examined so far (see also 
Joswiak &, Brownlee ( 200ll )). It is note- 
worthy that the IDPs show a 6.8 /im ab- 
sorption feature which is probably due to 



1 Based on observations with ISO, an ESA 
project with instruments funded by ESA Member 
States (especially the PI countries: France, Ger- 
many, the Netherlands and the United Kingdom) 
and with the participation of ISAS and NASA. 



the stretching vibration of carbonate an- 
ions (Sandford 1986). It has to be noted 
however, that there are doubts on this 
identification because the carbonate band 
at 11.4 /im was not detected in a study on 
five IDPs (Quirico et al. 2000). A band 
similar to the 6.8 fxm feature is supposedly 
seen in the spectra of the "Deep Impact" 
ejecta from comet 9P/Tempel 1 measured 
with the Spitzer IRS spectrometer (Lisse 
et al. 2006). In these spectra, also bands 
at 11-12 and 13-14 /xm have been assigned 
to carbonates. 

While the identificati on of the 1 1.3 ixm 
band with carbonates ( Bregman ( 19771 )) 
has been confuted by McCarthy, Forrest 
& Houck (1978) and subsequently been 
replaced by an identification with hydro- 
carbons (jCohen et al.l (119861 )). there has 
been continuous and recently even in- 
creased interest in carbonates as Stardust 
and planetary dust components. Carbon- 
ate bands at wavel engths larger than 2 5 fim 
came to the fore: Erard et al. (20001) and 



Lellouch. Encrenaz k, de Graauwl (l2000l ) 
suggested calcite and siderite (FeC0 3 ) as 
carriers of aerosol absorption bands of the 
Martian atmosphere inter alia at 26.5, 30- 
31 and 43.5 /im. 



Kemper et al.l (j2002al ) and lKemper et al 



(j2002bl ). by a careful analysis of the 2.4- 
200 /xm data of the Planetary Nebula NGC 
6302 taken with the Infrared Space Ob- 
servatory (ISO), were able to show that 
the previously unidentified '92' /im emis- 
sion band and the 62 /im component of 
the 60 /xm emission complex in this neb- 
ula can possibly be attributed to calcite 
and dolomite, respectively. According to 
their model fit, the carbonates belong to 
the cold (30-60 K) dust component of this 
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nebula. 



Chiavassa et al.l (120051 ) studied the 60 



180 /im ISO spectra of more than 30 pro- 
tostars and detected a solid state feature 
peaking between 90 and 110 /im in 17 of 
these spectra. Again, calcite has been sug- 
gested as the feature carrier, even though 
the accompanying characteristic resonance 
features of calcite could not be detected 
in the respective sources for different rea- 
sons (partly related to the sensitivity of 
ISO's detectors). For two star forming re- 
gions - the Carina Nebula and Sharpless 
171 - th is band had alre a dy be en discov- 
ered by lOnaka fc Okadal (l2003h but had 
been ascribed to onion structured carbon 
grains. The variation in the band position 
has been explained by Chiavassa et al. by 
varying amounts of other than calcium ions 
(e.g. Mg or Fe ions) in the carbonatic car- 
rier. 

The astronomical observations summa- 
rized so far indicate a need for more de- 
tailed investigations in two respects: (a) 
measurement of IR properties below K 
and (b) investigation of variations in the 
band profiles due to particle shape and size 
effects. Both is best achieved on the ba- 
sis of determining optical constants for low 
temperatures. Furthermore, this approach 
makes it possible to get rid of the possible 
influence of the matrix effect which may 
have influenced part of the previous com- 
parisons between laboratory and observa- 
tional data on carbonates. 

Given that it is unclear how carbonates 
can form in a circumstellar environment, 



Ferrarotti fe Gail! (120051 ) examined the for- 
mation of calcium carbonate in an oxygen- 
rich stellar wind (C/O < 1) by thermo- 
dynamic calculations, assuming the for- 



mation of calcite from Ca, H2O and CO. 
They came to the conclusion that calcite 
- if at all - would rather form in Asymp- 
totic Giant Branch (AGB) stars with low 
mass loss rates (approximately 1 x 10 -6 
M /y) than in AGB stars in their su- 
perwind phase (of which planetary nebu- 
lae like NGC 6302 are the remnants how- 
ever). But also in low-mass- loss-rate AGB 
stars, less than 1% of the total amount of 
du st would be carbo n ate du st. According 
to iFerrarotti fe Gail! (120051 ). it is mainly 



the relatively low condensation tempera- 
ture of carbonates which makes their con- 
densation in circumstellar environments so 
inefficient: for this temperature range to 
be reached in an expanding circumstellar 
shell, the dust must have a substantially 
larger distance from the central star (and, 
accordingly, a substantially reduced mate- 
rial density) than for the silicates. 

Toppani et al. (2005), however, sug- 
gest a possibly (more) efficient forma- 
tion mechanism of carbonates in the Post- 
AGB phase of stellar evolution. Based on 
laboratory experiments, they predict car- 
bonate formation by reactions between a 
comparatively hot silicate gas and a 300- 
500 K H 2 0-C02-rich gas under conditions 
far from local thermal equilibrium (LTE). 
They assume this non-LTE carbonate for- 
mation to take place when high-velocity 
winds of forming PNe interact with dense 
remnants of the slow-wind-AGB phase. 

A major goal of our study is to con- 
tribute to an assessment of the validity of 
the 'carbonate hypotheses' from the spec- 
troscopic point of view. 

The structure of the present paper is the 
following. In §2, we describe the deriva- 
tion of optical constants of CaCOs and 
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CaMg[C03]2 from reflectance measure- 
ments, both for room temperature and 
for cryogenic temperatures. In §3, ab- 
sorption efficiency spectra of small calcite 
and dolomite grains are presented. Fi- 
nally, in §4, we compare these small par- 
ticle spectra with infrared bands of as- 
tronomical objects which have been at- 
tributed to carbonates in recent studies. 
For the case of calcite, we present new 
constraints on the particle shapes and 
temperatures which are necessary to rec- 
oncile the '92' /xm b and assignments by 



tice structure of Mg and Ca carbonates 



sec 



Kemper et al.l (l2002al ) and lChiavassa et al. 



( 20051 ) with currently available laboratory 
data. 



2. Reflectance spectroscopy and de- 
rivation of the optical constants 

2.1. Crystallographic properties and 
sample preparation 

Calcite and dolomite both belong to the 
trigonal crystal class; hence, they are both 
uniaxial, and polarized reflectance mea- 
surements are required to determine their 
optical constantso The space group of 
calcite is R3c, while dolomite belongs to 
the space group R3 and has a reduced 
symmetry compared to its chemically sim- 
pler relative. While the calcium ions oc- 
cupy crystallographically equivalent sites 
in CaCOs, half of the places occupied by 
Ca ions in calcite are occupied by Mg 
ions in dolomite. The reduced symmetry 
of dolomite compared to calcite leads to 
the occurrence of an additional IR band 
(which is located close to 33 fxm) in its 
spectrum. For more details on the lat- 



. Hellwege et al.l f|1970r) an d references 
therein (e.g. IWyckhoffl (jl963 ) ) . 



2 Interestingly, the phenomenon of birefringence has 
been first observed in calcite crystals. 



In order to perform polarized reflectance 
spectroscopy, planes containing the respec- 
tive crystal's c-axes need to be found or to 
be produced. This is a difficult task in the 
case of ordinary CaC0 3 or CaMg[C0 3 ] 2 
rhombohedra, since none of their surfaces 
fulfill this condition. The c-axis runs di- 
agonally through the obtuse corners of 
cleaved rhombohedra. 

In the case of our dolomite rhombohe- 
dron - a 2x2x2 cm 3 crystal from Eugui, 
Navarra, Spain -, it was therefore neces- 
sary to cut it in order to produce a c- 
parallel surface, which has been polished 
subsequently. Polarized reflectance mea- 
surements on the polished surfaces showed 
no 'mixing' of the ordinary and extraordi- 
nary ray, which proved the cut to be paral- 
lel to the c-axis indeed. In the case of our 
calcite crystal - a 4 x 2 x 1 cm 3 fragment of a 
prism found at Sweetwater, Missouri, USA 
- no cut was necessary because this indi- 
vidual crystal had a relatively rare stature, 
its long edges being grown parallel to the 
symmetry axis. Hence, polishing one of 
the lateral surfaces of this crystal was suffi- 
cient to obtain a large measuring spot con- 
taining the c-axis. Indeed, polarized spec- 
troscopy on that spot did not show signif- 
icant traces of mixing of the two principal 
orientations of the electric field vector rel- 
ative to c (E||c, E_Lc; see also below, Fig. 

H). 

2.2. Reflectance spectroscopy 

The reflectance spectra of our two sam- 
ples have been measured at near-normal 
incidence using an FTIR spectrometer 
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Bruker 113v in the wavelength range 2- 
120 /im at a resolution of 2 cm -1 . Gold 
mirrors have been used as references. A 
turnable grid polarizer on polyethylene 
substrate positioned at the entrance of the 
sample compartment was used to select 
polarization directions parallel and normal 
to the plane of incidence. 

Fig. [1] shows reflectance spectra of cal- 
cite in the 5-120 /im range which comprises 
all IR bands of this mineral. In the case of 
the reflectance for the ordinary ray (E_Lc), 
an offset of 1.0 has been added to all re- 
flectance values. The different linestyles 
correspond to the different temperatures 
(10-300 K) as indicated. 

For the ordinary ray, five reflection 
bands can be discerned in the infrared 
spectrum of calcite, which are centered at 
6.8, 14.0, 30.7, 44.0 and 91.4 /an for room 
temperature. For the extraordinary ray, 
the respective reflectance band positions 
amount to 11.4, 29.6 and 91.4 /im (again 
for room temperature). 

In Fig. [21 the reflectance spectra of 
dolomite are plotted on the same scale as 
for calcite. Dolomite has a slightly more 
complex band pattern than calcite, which 
is a consequence of its comparatively re- 
duced lattice symmetry mentioned above. 
The reflectance maxima are located at 6.6, 
13.8, 25.5, 38.4 and 61.7/xm for the or- 
dinary ray (at T = 300 K), while for the 
extraordinary ray, the respective positions 
are 11.3, 25.4, 31.6, 56.7 and 60.0 /an. 

For the cryogenic measurements, the 
experimental setup was the following: 
the samples have been placed into a 
continuous-flow liquid-helium cryostat with 
contact-gas cooling (CryoVac KONTI Spec- 
tro B) equipped with polyethylene win- 



dows. It was inserted into the FTIR spec- 
trometer's sample compartment at the 
same position where the room tempera- 
ture reflection measurements had been per- 
formed. The cryostat windows did not al- 
low measurements at wavelengths shorter 
than 15 /mi. The cryostat possesses two 
sample mounts which can alternatingly be 
moved into the infrared beam, thus allow- 
ing subsequent measurements on the sam- 
ple and the reference mirror at each tem- 
perature. A third measurement with the 
sample being tilted to suppress the reflec- 
tion was necessary (again at each tempera- 
ture) in order to determine the reflectivity 
of the cryostat windows ("zero signal"), 
which had to be subtracted from both the 
sample and the reference spectra. The 
temperature was controlled by means of a 
CryoVac temperature controller using an 
ohmic heater and two silicon diode tem- 
perature sensors connected to the liquid- 
helium cooled walls of the contact gas 
chamber and to the sample mount, respec- 
tively. The temperatures have been stable 
to about ±0.5 K during the measurements. 

Spectra have been taken at 200 K, 
100 K, and 10 K during cooling down the 
samples. Two cooling cycles have been 
carried out for each sample with the two 
polarizer orientations. The wavelength- 
and temperature-dependent reflectances 
have been calculated by dividing the "zero- 
corrected" sample and reference spectra at 
each temperature. 

It is found that both for dolomite and 
calcite all bands in the above indicated 
wavelength range undergo a shift to shorter 
wavelengths, a bandwidth decrease and 
band strength increase with decreasing 
temperatures. These effects are especially 
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Fig. 1. — Polarized reflectance spectra of calcite for T = 300 K (solid line), 200 K (dashed), 
100 K (dash-dotted) and 10 K (dotted line). The spectra for E_Lc were shifted vertically by 
an amount of 1.0. 
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Fig. 2. — Polarized reflectance spectra of dolomite for the same series of temperatures as in 
the case of calcite. 
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significant for calcite's 44 /mi reflectance 
band and for dolomite's 38 /iin reflectance 
band. A more detailed discussion is given 
in §3. 

2.3. Lorentz oscillator fits and de- 
rived optical constants 

The polarized reflectance spectra of cal- 
cite and dolomite can be quite accurately 
reproduced by Lorentz os cillator fits, as 
has b een shown, e.g., by iHellwege et al. 
(Il970l ). Using the ansatz 



TO] 



where is the dielectric constant for 
frequencies large compared to all infrared 
wavelengths, Qj is the strength, TO., the 
transverse optical phonon frequency, 7^ the 
damping constant of the j-th oscillator, re- 
spectively. 

Figure [3] shows the optical constants of 
calcite which have been derived from the 
polarized reflectance spectra by means of 
Lorentz oscillator fits. The oscillator pa- 
rameters which were used for these fits - 
five oscillators for the ordinary ray and 
three for the extraordinary ray - are listed 
in Tab. U 

For dolomite, we only list its oscillator 
parameters - six oscillators for the ordi- 
nary ray and five for the extraordinary ray 
(see Tab. [2]) -, while the optical constants 
are not shown. However, both for calcite 
and dolomite, the optical constants can 
be retrieved from the electronic database 
http : / /www . astro . uni- j ena . de/Labo- 
ratory/OCDB , which also includes graphs. 



Apart from the different number of os- 
cillators, another important difference be- 
tween calcite's and dolomite's IR spectra 
is the absence of any phonon band for the 
latter at wavelengths larger than 67 /im (or 
wavenumbers smaller than 150 cm -1 ). In 
other words, it is noteworthy that calcite's 
IR bands extend to wavelengths roughly 
1.5 times larger than dolomite's 'reddest' 
band. 

For 300 K, our oscillator parameters are 
well comparab l e wit h those derived by 
Hellwege et al.l (11970l ) except for some mi- 
nor bands which are missing in their mea- 
surements (e.g. our '4:er for dolomite), 
probably due to a somewhat too small 
spectral resolution. 

3. Small particle spectra for differ- 
ent shape distributions 

From the optical constants shown and 
discussed in the previous section, the ab- 
sorption cross sections per particle volume 
- Cabs/V - have been calculated for grains 
of different shapes. Since the focus of this 
paper is on bands at large wavelengths, we 
choose the limit of particles being small 
compared to the wavelength (Rayleigh 
limit, size parameter x=27ra/A <C 1 and 
le 1 / 2 ! x < 1, a being the grain radius or 
largest semidiameter) . 

The following grain shapes have been 
taken into account: spheres (according 
to Mie theory) and continuous distribu- 
tions of ellipsoids (CDEs). As for the 
latter, two different modifications of it 
have been used: a CDE with equal proba- 
bility for all ellipsoid axis ratios ('mea n 
CDE', after iBohren fc Huffman! fll983h ) 
and a CDE with maximum probability 
for spherical grains ('weighted CDE', af- 
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Fig. 3. — Optical constants of calcite for room temperature, 200 K, 100 K and 10 K for the 
ordinary ray (E_Lc) (left) and for the extraordinary ray (E||c) (right). 
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Table 1: Lorentz oscillator parameters derived for calcite. 'OR' refers to the ordinary ray, 
'ER' to the extraordinary ray. Where no temperature is indicated as subscript, oscillator 
parameters have been derived from room temperature measurements. For the dielectric 
background, we chose the value e in f 5 oR = 2.6 in the case of OR and e^ER = 2.3 in the case 

of ER. All values are given i n units of cm -1 . 

j TOj Qj jj 



Ior 1407.4 1039.6 14.8 

2 OR 712.0 71.4 4.0 

3 O R300K 297.0 377.7 15.7 

3or2ook 302.5 385.4 12.0 

3oriook 306.0 387.0 8.87 

3oriok 307.1 387.6 7.34 

4 r3ook 222.8 206.1 11.0 

4 R2ook 226.5 211.2 8.35 

4oriook 228.8 219.1 6.14 

4 riok 229.7 219.7 4.54 

5 O R300K 104.7 148.1 6.00 

5 R2ook 103.2 165.7 5.35 

5oriook 103.8 158.0 3.57 

5qr,iok 104.2 167.6 2.45 

\ E r 872.0 250.8 2.09 

2er,30M 306.7 368.0 16.9 

2er,2ook 312.6 372.2 11.6 

2 E r,iook 317.6 377.6 9.19 

2er,wk 318.9 378.7 7.30 

3 E r,30ok 94.3 200.6 7.10 

3 E r,20ok 93.3 203.2 5.91 

3er,wok 94.6 203.6 4.37 

3eraok 95.0 202.4 2.86 
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Table 2: The set of Lorentz oscillator parameters derived for dolomite by fitting the room 
temperature reflectance spectra. All designations have the same meaning as in the previous 
table, but here, we chose e in f 5 oR = 2.7 and e infjER = 2.4. 
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TOj 


Qj 


7j 


loR 


1453.8 


399.2 


23.8 
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1425.5 


980.0 


15.0 


3or 
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98.8 


6.14 


4oR,300K 
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4oR,200K 
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2 ER 300 K 


364.8 


335.8 


20.4 


2 ER,200K 
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3eR,300K 


316.8 


179.5 
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319.4 
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8.24 


3er,wok 
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ter lOssenkopf. Henning &: Mathia (119921 ) ) . 



3.1. Temperature effects for spheri- 
cal grains 

Figures H] and [5] show the volume nor- 
malized absorption cross sections C a & s /V 
calculated for spherical calcite and dolomite 
grains, respectively, depending on temper- 
ature. The figures demonstrate that the 
effect of the cooling on the band posi- 
tions and band widths is present also in 
the small particle spectra. Similarly to the 
reflectance spectra, this effect is especially 
pronounced for the case of the ~40 /im 
band of calcite, which shifts from 42.6 /zm 
for 300 K to 41.4 /im for 10 K, while at the 
same time, the C a fe s /V peak value increases 
by a factor of 2.4. For the calcite bands 
around 30 and 80 /im, the temperature de- 
pendence is not equally strong, but still 
significant. 

In the case of CaMg [003)2, the tempera- 
ture-dependent evolution of the spectra is 
similar apart from the effect that there is 
only an insignificant change between 100 K 
and 10 K (cf. Fig. [5]). For the additional 
band in the dolomite spectrum (compared 
to calcite) around 31 /im, the temperature 
dependence is rather moderate compared 
to the other bands. In the wavelength 
range A < 15 /zm, the dolomite and the 
calcite spectra resemble each other very 
closely. 

3.2. Grain shape effects at low tem- 
peratures 

Due to the high oscillator strengths 
characterizing most of the infrared bands 
of both calcite and dolomite, the grain 
shape dependence of the band profiles is 
very strong. For the FIR bands of cal- 



cite, this is evident from the comparison 
of Fig. H] with Fig. [6] as well as from the 
peak wavelengths given in Tab. [3j Cal- 
cite's 30 /im band is shifted to the 'red' 
by about 2 /im for a mean ODE compared 
to the case of spherical grains, and at the 
same time its width is increased by factor 
of four. The 80-90 /im double-band of cal- 
cite degenerates into a broad single bump 
centered at 90.1 /zm with two wings for the 
mean ODE. For the weighted ODE case, 
the maximum of the absorption efficiency 
is located at 86.3 /im for a grain temper- 
ature of 300 K. For the 40 /im band, the 
broadening due to the shape distribution 
is not equally strong as for the 80-90 /im 
band, and the shift of the band position to 
longer wavelengths is limited to 1 /im for 
a given temperature (compared to more 
than 5 /im for the 80-90 /im band). 

Comparing the band properties for dif- 
ferent temperatures within a given grain 
model, we find that the shift of the peak 
position is very prominent in the case of 
the 40 /im feature, even for the grain shape 
distributions. It clearly dominates over 
shape effects for this band, whereas this is 
not true for the 30 and 80-90 /im features. 
In the case of the latter both the shape 
effects are stronger due to the higher os- 
cillator strengths and the temperature ef- 
fects are intrinsically smaller as has been 
discussed in the previous section. 

3.3. Comparison to PE powder spec- 
tra 

Since powder spectra of calcite and 
dolomite - measured at room tempera- 
ture - have been utilized by both Kemper 
et al. (2002a) and Chiavassa et al. (2005) 
in order to argue for the presence of car- 
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Fig. 4. — Volume normalized absorption cross sections of spherical calcite grains for room 
temperature, 200 K, 100 K and 10 K in the wavelength range 5-110 //m. 
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Table 3: Positions of the maxima of C a b s /V for calcite grains at different temperatures and 
for different shape distributions 





20 /iin complex 


40 jum 


80-90 iim 


Peak position 


H 


H 


H 


sphere 








10 K 


27.05, 27.84 


41.40 


79.62, 83.77 


300 K 


28.58 


42.64 


80.56, 84.71 


weighted CDE 








10 K 


28.18 


41.75 


84.91 


300 K 


29.13 


43.09 


86.25 


mean CDE 








10 K 


30.02 


42.54 


90.31 


300 K 


30.66 


43.68 


90.14 




300K, mean CDE 
10K, mean CDE 
300K, weighted CDE 
10K, weighted CDE 



A [ftm] 

Fig. 6. — Cabs/V profiles for a distribu- 
tion of ellipsoidally shaped calcite grains 
(see text for details). Only the wavelength 
range 20-120 /im and only the two tem- 
perature 'extremes' T=10K and T=300K 
have been depicted here. 



bonate dust in astronomical objects, it is 
important to discuss the relation between 
such powder spectra and the small particle 
spectra based on our optical constants. By 
'powder spectra' or 'PE spectra', we un- 
derstand the transmission spectra of sub- 
micron sized ground powder embedded in 
polyethylene (PE), transformed into ab- 
sorption cross sections. 

Powder spectra generally depend on 
various circumstances such as particle size 
and shape, orientation, agglomeration de- 
gree, and the refractive index of the sur- 



rounding medium (e.g. lBohren fc Huffman 
fll983f n. 

Fig. [7] shows both the powder spectrum 
published by Kemper et al. (2002a) - mea- 
sured on a sample of unsedimented calcite 
grains (dash-dotted-line) - and of a sam- 
ple of calcite grains which had previously 
been subjected to sedimentation in a liquid 
(acetone), in order to remove large grains 
(solid line). The peak of the latter spec- 
trum is shifted by about 2 fim to shorter 
wavelengths with respect to the data pub- 
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lished by Kemper et al. (2002a), although 
both have been measured at room tem- 
perature. This is due to the presence 
of large grains (several /im in diameter) in 
the unsedimented powder. The spectrum 
of the sedimented powder is in good agree- 
ment with the result of the calculation for a 
mean CDE with PE environment (see Fig. 
Cj), indicating the presence of a consider- 
able fraction of grains with shapes far away 
from spherical symmetry. 

Conversely, we find that the mean CDE 
is a good representation of the shape distri- 
bution of a real calcite powder with grain 
sizes within the Rayleigh limit. This jus- 
tifies the use of this shape distribution 
for the comparison with observed data at 
large wavelengths. In the following (see 
Sect. 4.1) we will apply this model to sim- 
ulate dust spectra at low temperatures. 
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Fig. 7. — Comparison of calcite's FIR 
bands measured by PE powder transmis- 
sion spectroscopy (solid and dashed-dotted 
lines; see text for details) with the absorp- 
tion efficiency spectrum calculated from 
the optical constants for a mean continu- 
ous distribution of ellipsoids (dashed) and 
for a weighted CDE (dotted), both in PE. 



4. Comparison with astronomical 
observations 

The laboratory spectra presented so far 
are especially relevant for a comparison 
with FIR spectra of cold dust in objects 
such as the planetary nebula NGC 6302. 
Hence, we shall first give a brief review of 
some key properties of the dust spectrum 
of this object. 

Barlow (1997) observed features 65 and 
69 /im and another possible broad band 
feature in the range of 88-98 /mi. While 
the 65 /im feature was attributed to crys- 
talline ice (see also Waters et al. 1996), the 
69 /im feature was suspected to be due to 
crystalline forsterite and the question for 
the carrier of the 88-98 //m band remained 
open. 

Molster et al. (2001) identified the 
blend of diopside and 
crystalline water ice with enstatite. They 
also confirmed the reality of the broad 
feature around 90 /im and suspected that 
there should be a very cold dust compo- 
nent present in the nebula. 

Kemper et al. (2002a,b) first assigned 
the ~65 fim feature to cold dolomite (with 
additional contributions by diopside and 
water ice) and the broad emission band 
around 90 /im to cold calcite dust. As al- 
ready mentioned above, the temperature 
of the carbonate dust was assumed to be 
in the 30-60 K range. A mass fraction of 
less than 0.3% (for calcite as well as for 
dolomite) has been derived by these au- 
thors. 

We retrieved an ISO-SWS and an ISO- 
LWS spectrum of NGC 6302 from the ISO 
archive and reduced it by means of the 
OLP version 10.0. 
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From the composite spectrum (ranging 
from 2.4-200 /xm), we subtracted a com- 
bination of Planck functions for tempera- 
tures of 30, 55 and 95 K. The remaining 
'residual dust emission' is shown in Fig. 
[8] for the 35-105 fim range; the individual 
broad peaks of the gray line result from 
various ice and dust species such as H 2 0, 
forsterite, enstatite, diopside, and - possi- 
bly - carbonates (see Kemper et al. 2002a, 
Tab. 2). 



NGC 6302 (cent, subtracted) 

calcite, CDE, 30K (abs.eff. x Plonck function) 




A [p] 



Fig. 8. — Comparison between the residual 
dust emission of NGC 6302 and synthetic 
dust emission spectra (for 30 K and 60 K) 
of calcite and dolomite in the 35-105 fim 
range. 

Our aim has been to examine whether 
the assignments of the 60-65 /im and '92' fim 
bands to dolomite and calcite is compati- 
ble with the optical constants and powder 
spectra presented in Sect. 3, and, if so, for 
which temperatures and particle shapes 
this is the case. For answering this ques- 
tions, the detectability of bands at wave- 
lengths between 35 and 45 fim depending 
on the dust temperature also plays a key 
role. 



4.1. The '92' fim band in NGC 6302 

A detailed view of the '92' /im band 
is given in Fig. [9J It can be seen from 
this figure (when compared to Fig. [6]) that 
only for a mean CDE particle shape dis- 
tribution of calcite grains, the bandwidth 
is as large as in the case of the NGC 
6302 spectrum. For a weighted CDE, the 
bandwidth is smaller by 45% (at the rel- 
evant cryogenic temperatures). The best 
fit is achieved by the C a b s /V function for a 
mean CDE multiplied by a Planck func- 
tion for a temperature of 30 K. A dust 
temperature of 60 K is hardly compatible 
with the observations, both with respect 
to the '92' /im band profile and with re- 
spect to the emergence of a ~42 /im fea- 
ture (see dotted lines in Figs. M and [9]). 
The calcite data derived from transmission 
spectroscopy of sedimented powder (dash- 
dotted line in Fig. [9]) roughly reproduce 
the observed band profile as well. It is 
noteworthy that the '92' /im band profile 
of NGC 6302 obtained by Kemper et al. 
(2002a) is better fitted by the PE spec- 
trum of an unsedimented calcite powder 
sample which is broader and peaks at a 
larger wavelength. This may be due to the 
use of different ISO LWS spectra (a mean 
of seven vs. a single) or due to some uncer- 
tainty in the continuum subtraction. 

4.2. The 60-65 fim emission com- 
plex in NGC 6302 

In Fig. [10] a close-up view of the 60 /im 
emission complex is shown. Although the 
bandwidth of a mean CDE particle shape 
distribution can reproduce the observed 
bandwidth of the 60 /im feature, neither 
for a mean CDE nor for a weighted CDE - 
both at a temperature of 45 K - the posi- 
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Fig. 9. — Comparison between the residual 
dust emission of NGC 6302 and synthetic 
dust emission spectra of ensembles of ellip- 
soidally shaped calcite grains for 30 K and 
60 K in the 75-105 /mi range. 



Fig. 10. — Comparison between the resid- 
ual dust emission of NGC 6302 and syn- 
thetic dust emission spectra of ensembles 
of ellipsoidally shaped dolomite grains for 
30 K and 50 K in the 49-81 /mi range. 



tions coincide with the observed one. The 
difference between the band positions re- 
sulting from laboratory measurements and 
the observed band position (~ 64 /mi) is 
too large as to be compensated by tem- 
perature effects - see also Fig. El where a 
comparison with dolomite dust at T=30 K 
is made. Dust temperatures lower than 
40 K would allow to suppress the dolomite 
emission band at ~37/im (see again Fig. 
[S]) , but the position of the 60 /an feature 
would still not match. Therefore this fea- 
ture cannot be fit only by dolomite, while it 
is not excluded that dolomite contributes 
to the 'blue wing' of the 60 /mi emission 
complex. Our study confirms that the 60- 
65 /im band must have other major car- 
riers, e.g. diopside as suggested by Kem- 
per et al. (2002a). A recent publication by 
Chihara et al. (2007) shows that a mixture 
of diopside and akermanite (Ca 2 MgSi 2 07) 
at a temperature of 30 K provides a good 
match of the 60-65 /xm feature. 



4.3. The 90-110 /im band in proto- 
stars 

Chiavassa et al. (2005) studied a sam- 
ple of 32 low and intermediate mass pro- 
tostars (mostly Class and Herbig Ae/Be 
objects). Out of this sample, 17 sources 
show a broad band feature between 90 and 
110 /im. The dust temperature in these ob- 
jects is found by blackbody fits to the con- 
tinuum to be below 50 K (except for emis- 
sion from the innermost regions). Calcite 
is suggested to be a possible carrier of this 
feature, which seems to be supported by 
the laboratory measurements of Kemper et 
al. (2002a) that show a peak at 93 /im and 
a width (FWHM) of 16 /mi. 

Our data show that for PE powder mea- 
surements on small grains and also for cal- 
culated spectra assuming a mean CDE, the 
maximum of the '92 '/mi calcite band is lo- 
cated at 90/mi. For other particle shape 
distributions, it is narrower and located at 
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much smaller wavelengths, especially for 
the temperature range relevant for PNe 
and protostars. The transverse optical 
resonance wavelength of mode 5or,wk of 
about 96 /im (Tab. 1) can be considered as 
the maximum peak pos i tion f or this band 
(cf. iBohren fc Huffman! ril983h h The fea- 
ture produced by the 3er,wk mode with a 
TO resonance wavelength of about 105 fim 
is much broader and therefore compara- 
bly weak and cannot dominate the cal- 
cite grain emission even if alignment effects 
were present. 

Hence we conclude that the bands ob- 
served for the Herbig Ae/Be objects may 
in principle be compatible with emission 
from calcite grains, if sufficiently extreme 
shape distributions (much more extreme 
than the mean CDE) and/or large grains 
are assumed (cf. Fig. 7 in Chiavassa et al. 
2005). For the Class objects, in the con- 
trary, we confirm that the feature charac- 
teristics are incompatible with laboratory 
data for calcite. Whether the incorpora- 
tion of other ions into a carbonate can shift 
the band to still larger wavelengths has to 
be left open at this point. 

5. Conclusions 

While previous studies on carbonate 
dust in astronomical environments were 
largely based on powder measurements 
performed at room temperature, we have 
on the one hand derived optical constants 
of calcite and dolomite and on the other 
hand carried out measurements at cryo- 
genic temperatures (200, 100 and 10 K). 
From these optical constants, small parti- 
cle spectra (absorption efficiencies as func- 
tions of the wavelength and temperature) 
were calculated for different grain shapes. 



Since the optical constants of calcite and 
dolomite reach peak values of 5-10 (corre- 
sponding to maxima of the imaginary part 
of the dielectric function of up to 150), 
their absorption efficiencies are extremely 
sensitive to particle shape effects. For the 
'92' /im band of calcite, e.g., the negativity 
range of the real part of the dielectric func- 
tion is as large as 73-96 /im for the ordinary 
ray (and even 75-105 /im for the extraor- 
dinary ray). The former values approxi- 
mately indicate how large the range is in 
which calcite's '92' /im band can shift due 
to particle shape effects, since the ordinary 
ray is weighted more strongly than the ex- 
traordinary ray. 

At cryogenic temperatures, the FIR 
bands of calcite and dolomite increase 
their intensity, get narrower and shift to 
shorter wavelengths. The temperature- 
related shift in position is largest for cal- 
cite's 40 /im band: in this case, it reaches 
0.5 /im per 100 K. 

As for the '92' /im feature detected by 
Kemper et al. (2002a) in ISO-LWS spectra 
of NGC 6302, this feature can only be as- 
signed to CaC0 3 grains if the assumption 
of a mean CDE - with equal probability of 
all grain shapes - is made. 

Furthermore, dust temperatures have to 
be well below 40 K in order to suppress the 
30 /im and 42 /im bands of calcite which 
accompany the '92' /im band but are not 
seen in the spectra of NGC 6302. 
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of NGC 6302, reduced with the ISO spectral 
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